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ABSTRACT 

Silicon monoxide maser emission has been detected in the circumstellar envelopes of 
many evolved stars in various vibrationally-excited rotational transitions. It is consid- 
ered a good tracer of the wind dynamics close to the photosphere of the star. We have 
investigated the polarization morphology in the circumstellar envelope of an AGB star, 
R Cas. We mapped the linear and circular polarization of SiO masers in the v = 1, 
J = 1—0 transition. The linear polarization is typically a few tens of percent while the 
circular polarization is a few percent. The fractional polarization tends to be higher 
for emission of lower total intensity. We found that, in some isolated features the frac- 
tional linear polarization appears to exceed 100%. We found the Faraday rotation is 
not negligible but is ~ 15°, which could produce small scale structure in polarized 
emission whilst total intensity is smoother and partly resolved out. The polarization 
angles vary considerably from feature to feature but there is a tendency to favour the 
directions parallel or perpendicular to the radial direction with respect to the star. In 
some features, the polarization angle abruptly flips 90°. We f ound that our data are 
in the regime gfl ^> R 3> T, which indicates that the model of iGoldreich et ail (|l973l) 
can be applied and the polarization angle flip is caused when the magnetic field is at 
close to 55° to the line of sight. The polarization angle configuration is consistent with 
a radial magnetic field although other configurations are not excluded. 
Key words: maser - polarization - magnetic field - star: AGB - star: late-type star: 
individual: R Cas . 



1 INTRODUCTION 

The circumstellar envelope (CSE) in an asymptotic giant 
branch star (AGB) is a very active region. Any AGB star 
will lose mass in the form of a slow wind at a rate that 
will significantly affect the mass of the star as well as en- 
riching the interstellar medium with nuclear processed ma- 
terial. This produces a circumstellar envelope of escaping 
dust and gas particles. In the zone just above the photo- 
sphere, out to 57?* (stellar radius), the extended atmosphere 
experiences periodic shocks, so each period is divided into 
intervals of mass outflow and infall. In this zone, the cycle- 
averaged temperature drops as the wind flows away from 
the star, f rom about 3100 K at the photosphere to « 750 K 
at 5.57?+ (|Grav et al.ll2009h : cooling takes place by line ra- 
diation from various molecules, especially H2O. SiO masers 
provide sub-milliarcsec images of this region. 

Some molecular species are formed in an equilibrium 
process deep in the atmosphere and are destroyed in the 
outer parts of the outflow by interstellar ultraviolet radiation 
(H2, CO, H2O). Other species, for example SiO, are depleted 



due to conden sation on dust particles at a few stellar radii 
l|Habing|ll99rJ) . 

The SiO maser zone is thought to be shaped by a com- 
bination of shocks from stellar pulsations, gravity and possi- 
bly magnetic effects, with radiation pressure on dust becom- 
ing possible towards the outer region. SiO masers are seen 
in rings, wi th faint or no emission along the line of sigh t 
to the star (|Diamond fc Kemball|[2u03l : ICotton et.al.ll2006h . 
This is due to tangential beaming from an accelerating, ap- 
proximately spherical wind, where masers have the deepest 
amplification paths if they are at similar velocities and in a 
similar plane of the sky to the star. 

Studying the magnetic field properties is done by in- 
vestigating the polarization morphology of the maser emis- 
sion. The maser emission is significantly linearly polar- 
ized but circular polarization i s weaker, as expec ted for 
a non-paramagnetic molecule (|Herpin et al.l 120061 ). Frac- 
tional circular polarization of SiO masers in late-type star s 
from single-dish observations is m c < 0.5% (|Habinelll996h 
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while the degree of linear polarization is me ~ 15 — 30% 
jTroland et al.lll979l ). 

There are two explanations suggested for the linear po- 
larization of SiO masers. Polarization can be produce d by 
the magnetic field (|Goldreich et al.lll973l ; lEhtzurll 19921 ). Al- 
ternatively, in the absence of a magnetic field, anisotropic 
pumping can produce s trongly polarized maser emission 
\ Western fcWatsonl 1 19831 ). This mechanism has been pro- 
posed as the cause of the tangentially polarized maser emis- 
sion seen in the map of SiO masers observed by the VLBA 
jDesmurs et.alJl200Ch . 

If it is magnetic in origin, the linear polarization posi- 
tion angle provides information about the structure of the 
magnetic field in circumstellar envelopes. The theoretical re- 
lationship between the polarization position angle and the 
projected magnetic field direction depends on 8f, the angle 
between the magnetic field direction and the maser prop- 
agation direction (our line of sight). When dp < 55°, the 
linear polarization vectors are parallel to the magnetic field 
lines and when 8f >55° the vectors are perpendicular to the 
field lines. When 8f ~55° the linear polarization vectors can 
flip ^ within a single feature. At this value of 8f, the Van 
Vleck angle of 55°, the masing region reaches maximum in- 
tensity and the fractional linear pola rization approaches zero 
(|Elitzuilll992l ; iGoldreich et alJll973h . 

R Cassiopeia is an oxygen-rich AGB star which is classi- 
fied as an M-type Mira-variable. Its optical brightness varies 
from magnitude +4.7 to + 13.5 with a period of 430 days and 
its mass is about 1.2 Mr^ . IVlemmings et alj (|2003l ) used as- 
trometric VLBI to measure a distance of 176^4g pc with a 
proper motion of (85.5 ±0.8, 17.5 ±0.7) mas yr" 1 in R. A. 
and Dec, respectively. Various estimates of the stellar ve- 
locity v-,, appear in the literature, we have adopted the value 
of 24±2 km s _1 . 

We monitored the 43 GHz SiO masers of R Cas for 
almost two stellar cycles. Th e total intensity results were 
discussed in I Assaf et al.l(|201ll ). Here, we describe the polar- 
ization observation in section [2] We present the polarization 
results in section [31 and analyse these results in section 3] 
drawing conclusions in section [5] 



2 OBSERVATIONS AND DATA REDUCTION 

The SiO masers around R Cas were observed as part of a 
more extensive programme of VLBA0 monitoring of other 
stars. In total we have 23 epochs of data for R Cas. Data 
were recorded at each VLBA antenna in dual-circular po- 
larization in two 4 MHz windows, each digitally sampled at 
the full Nyquist rate of 8 Mbps in 1-bit quantization. The 
lower spectral window was centred at a fixed frequency cor- 
responding to the ii = 1, J = 1 — SiO transition, at an 
assumed rest frequency of 43.12207 GHz and a systemic ve- 
locity Vlsr = +24 km s _1 . R Cas was observed for three 
45-minutes periods evenly spread over the 8 hour duration 
of the run. Adjacent to each R Cas observation, 5 minutes 
was spent observing the continuum calibrator 0359+509 at 

The VLBA (Very Long Baseline Array) is operated by the the 
National Radio Astronomy Observatory, a facility of the National 
Science Foundation operated under a cooperative agreement by 
Associated Universities, Inc. 



the same frequency as R Cas. The data were correlated at 
the VLBA correlator in Socorro, NM. The correlator ac- 
cumulation interval was set to 2.88 seconds. All polariza- 
tion correlation products (RR, RL, LR, LL) were formed. 
This configuration produced auto- and cross-power spectra 
in each 4 MHz baseband with a nominal velocity spacing of 
~ 0.2 km.s" 1 . 

We reduced the data using the standard approach 
to VLBI spectroscopy within the NRAO AIPS package 
(http://www.aips.nrao.edu/cook.html). We processed the 
visibility data using the the semi-automated spectral- 
line polarization calibration pipeline which was origi- 
nally w ritten by Diam o nd in 1998, based on the formal- 
ism in iKemball et al. I (1 19951) and developed further by 
Kemball fc Diamond! (1 19971 ) (for more detail see l Assaf et al.l 
20111) . 

After the pipeline processing, we noticed that there were 
systematic offsets in flux density of a few to 15% between 
LL and RR spectra for each epoch. We assumed that the 
net circular polarization of SiO across the band is zero and 
the total intensity calibration is correct. This enabled us to 
calculate corrections to align the flux scales for LL and RR 
and obtain realistic values of Stokes V, albeit with large 
uncertainties. 

The absolute Electric Vector Position Angle (evpa) for 
any linearly polarized emission is unknown because there is 
no instrumental measurement of absolute R-L phase differ- 
ence in the reference antenna of the VLBA and in the VLA, 
but the VLA can provide absolute astronomical calibration 
of evpa relative to a small subset of primary astronomical 
calibrators, given the lower spatial resolution of the array. 
We used the VLA to transfer absolute polarization calibra- 
tion to VLBA data via the compact secondary polarization 
calibrator (0359+509) which was observed by both arrays 
oil 19 th Dec 1999. The same reference antenna kp was used 
in the reduction of all epochs of VLBA data and the po- 
larization angle stability is likely to be better than 35° (C. 
Walker private communication). Hence, we were able to es- 
timate the polarization angle of 0359+509 as 30° ± 10°. 

We fitted 2-D Gaussian components to each patch of to- 
tal intensity (Stokes /) maser emission brighter than 5cr rms , 
to measure the p osition and flux density, as described in 
lAssaf et al.1 (|201ll ). We measured the Q, U and V flux den- 
sities at the position of the peak of the I component, as 
described in lKemball et all (|2009l ) . We use the peak (rather 
than integrated) intensity measurements of Stokes I for com- 
parison with the other Stokes parameter measurements. We 
estimated the uncertainties based on (beam si ze) / (signal-to- 
noise ratio), as appropriate for a sparse array (jCondon et.al.1 
1 19981 , iRichards et aHll999h . 

We formed the linearly polarized intensity 

P = (Q 2 + U 2 f /2 (1) 
the fractional circular polarization 

V 

m c = j (2) 
the fractional linear polarization 

mt = Y (3) 
and the electric vector polarization angle evpa 
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EVPA = 0.5arctan^- (4) 

(e.g. as defined bv lHeileslliooi or lElitzunll992T ). 

We measured the angular FWHM (s) of the total inten- 
sity components by deconvolving the restoring beam from 
the measured size. This was used to estimate the bright- 
ness temperature Tb from the integrated intensity. Series 
of components form features and we estimated the average 
positions, angular size L, angular FWHM d,F, peak flux den- 
sities and other properties of each feature. 

The largest angular separation between components 
making up a feature, L, measures the actual angular size 
of the emission detected. The separation of the components 
with flux density closest to half the peak, cfc , represents the 
beame d size of the feature, as explained in [Richards et alj 
|201ll ). with reference to Elitzur (1992). The beaming angle 
dO, = (^-) 2 . We note that cLf and L may be underestimated, 
as figs Bl-3 in lAssaf et al.l l|2011f ) show that a significant 
amount of total intensity flux is resolved out. There is no 
clear, systematic difference between the effect on brighter 
or fainter emission, so we suggest that the main effect is to 
increase the uncertainty in dfi, but it could be somewhat 
overestimated if the fainter emission is more diffuse, leading 
to L being more underestimated. 

An improved method of finding the centre of emission 
has le d to some minor changes with respect to lAssaf et alJ 
(2011), which do not affect those results significantly. Ta- 
ble flj shows the updated values of the average shell radius 
and the shell thickness, and the difference with the old val- 
ues. The differen ces are close to the uncertainties (given in 
lAssaf et al.ll201ll) . which are higher for the last few epochs 
due to the fainter emission and large gaps in the maser shell. 



3 RESULTS 

We made Stokes I, Q, U and V image cubes of R Cas, at each 
epoch, at spatial and spectral resolutions of approximately 
(40 x 20) /ias 2 , 0.2 km s" 1 . 

The polarization detection threshold for individual com- 
ponents is 5cr rms and me > 5% or m c > 15%; lower thresh- 
olds are possible when averaging over larger spectral or spa- 
tial regions (including correction for Ricean bias). Each fea- 
ture is made up of many components, and every feature con- 
tain more than 0.5 Jy/beam summed total intensity. Quiet 
channels have cr rma 0.03 — 0.04 Jy beam -1 , so any linear po- 
larization with me 20 percent will also have a signal to 
noise ratio 5. 



3.1 Circular Polarization 

The mean degree of circular polarization for each epoch was 
estimated using Eq. [2] 

We found the fractional circular polarization is m c ~ 
0.4 — > 6%. W e attempted to fit the expression given by 
lElitzurl (|l996h to the first derivative of Stokes V. This failed 
due to the weakness of Stokes V intensity and the large chan- 
nel separation which is greater than the Zeeman splitting. 



Epoch 


New (R) 


error 


dr 


Abs(diff) 


BD62A 


25.93 


0.67 


8.72 


1.16 


BD62B 


25.49 


0.05 


9.05 


1.59 


BD62C 


24.72 


1.77 


8.21 


1.33 


BD62D 


24.24 


1.22 


8.60 


1.81 


BD62E 


24.86 


1.59 


7.87 


2.24 


BD62F 


25.45 


1.98 


5.80 


2.98 


BD62G 


25.96 


2.11 


6.20 


1.45 


BD62H 


30.34 


0.30 


7.92 


2.25 


BD62I 


27.50 


0.53 


5.38 


0.99 


BD69A 


26.77 


1.64 


7.52 


0.52 


BD69B 


27.97 


0.37 


6.19 


0.26 


BD69C 


27.55 


1.50 


5.68 


0.45 


BD69D 


28.09 


0.16 


5.51 


0.11 


BD69E 


28.02 


0.50 


6.07 


0.52 


BD69F 


28.32 


1.71 


6.11 


0.17 


BD69G 


24.24 


2.50 


8.23 


4.28 


BD69H 


29.99 


0.26 


9.62 


0.83 


BD69I 


24.90 


0.17 


7.68 


2.25 


BD69J 


24.60 


0.63 


14.05 


2.23 


BD69K 


21.41 


3.11 


13.51 


3.70 


BD69L 


18.17 


1.22 


4.50 


5.79 


BD69M 


17.27 


3.83 


7.82 


8.60 


BD69N 


21.97 


1.68 


11.21 


4.85 



Table 1. The new shell radius (New i?),the error in R, the shell 
thickness (dr) and the absolute difference between the old and 
the new values, all in mas. 



3.2 Linear Polarization 

Figs. [T] [2] and [3] show the polarization morphology of SiO 
masers in R Cas for 23 epochs. The linear polarization vec- 
tors are superimposed on the map of the total intensity. The 
orientations of the vectors define the electric field plane and 
the length of the vectors are proportional to the linearly 
polarized intensity P. Fig. U shows an enlargement of the 
components making up a single feature. The magenta cross 
marks the error-weighted centroid of the feature. The thin 
magenta circle encloses half the feature flux, i.e. its diame- 
ter is the feature FWHM, c?f and the diameter of the thicker 
magenta circle represents the total feature extent L. 

We found that SiO maser emission is significantly lin- 
early polarized. We calculated the fractional linear polar- 
ization using Eq. [3] We found that the percentage of lin- 
ear polarization (averaged over all features per epoch) is 
rri£ ~ 11 — > 58% but exceeds 100% in some isolated fea- 
tures, discussed in Section f4. 21 

We investigated the relationship between the evpa and 
the radial direction with respect to the star, defined by the 
position angle in the plane of the sky, 8. Fig.[S]shows the pro- 
portion of the polarized emission within bins of (evpa— 8) 
in the whole SiO shell and in the inner shell (within the ra- 
dius enclosing 25% of the total maser flux at each epoch). 
The thickness of the line is proportional to the logarithm 
of the total linearly polarized flux. During the first cycle 
(<f> — 0.158 to <f> = 1.158, 15 epochs), we found that the po- 
larized flux in the whole shell is dominated by emission with 
evpa parallel to 8 at 10 epochs. The most popular evpa ori- 
entation with respect to the position angle 8 is perpendicular 
for the other 5 epochs, dominating in 3 of these. However, 
in the inner shell, only half the epochs have a high propor- 
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Figure 2. Same as Fi^T]but for Epochs 8-15 
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Figure 3. Same as FijJT]but for Epochs 15-23 
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tion of emission with parallel (6 epochs) or perpendicular (2 
epochs) EVPA; it is at intermediate angles for the remaining 

7 epochs of the first cycle. In the second cycle, 7 out of the 

8 epochs (from (f> = 1.234 to (f> = 1.783), are dominated by 
emission with intermediate evpa in the whole shell but this 
comes mainly from a single feature. At cj> = 1.662 the EVPA 
is predominantly parallel. The inner shell has more emission 
with a parallel evpa at epochs cj> — 1.501, 1.662, otherwise 
behaving similarly to the whole shell. 

The polarization structure during the first stellar cycle 
can be summarized as a bimodal distribution. Most of the 
linear polarization vectors are either radial (parallel to the 
position angle of the location of the emission in the projected 
shell) or tangential (perpendicular). However, the polariza- 
tion in the inner part of the shell is somewhat less ordered. 
The later parts of the second stellar cycle do not show any 
clear pattern, but the emission generally was noisier with 
fewer significantly polarized features. 




R.A. offset (mas) 



4 DISCUSSION 

We noted in Section [T] that there are two competing mod- 
els put forward to explain SiO polarization. The adopted 
model for the transfer of the polarized maser emission in a 
non-paramagnetic molecular transition plays an important 
role in the int erpretation of the polarization measurements 
|Elitzurlll996h . 

First, we made an independent estimate of the mag- 
netic field strength to test the feasibility of the magnetic 
model. The bulk kinetic energy density (i?Buik = ^P v ' 2 ) is 
^ 2.1 x 10~ 3 J m -3 , where p is the volume dens ity of SiO, 
taken as 1.42xl0" 13 gm cm -3 l|Grav et al.ll2009f). and v is 
the e xpansion velocity which is up to 5.5 km s _1 ()Assaf et alJ 
l201l| y The thermal energy density -Exhcrmai = ^nu 2 kT is 
0.832 x 10~ 3 J m~ 3 , where riH 2 is the molecular hydrogen 
density, k is Boltzm ann's constant an d T is the effective tem- 
perature, 1330 K l|Grav et al.ll2009l ). A first estimate of B 
can be found by comparing the energy densities. The mag- 
netic energy density Eb is given by 



2/Xo 



4 x 10" 



J m 



(5) 



By equating these en ergy densities , we fo und the first esti- 
mate of B ~ 0.725 G. iHerpin et al.l (|2006h measured a mag- 
netic field strength in the range 0.9 — ¥ 2.8 G from single dish 
observations of R Cas at 86 GHz, slightly higher than our 
estimate. 



4.1 Application of Zeeman interpretation 

We can decide which interpretation to apply by comparing 
three parameters: the stimulated rate emission R, the Zee- 
man splitting rate gQ, and the decay rate V. 

We calculated the maser beaming angle d£l and the 
maser brightness temperature Tb from our data. We found 
the error-weighted beaming angle (Table [2} averaged over 
23 epochs is dQ ~ 0.047 sr and the mean brightness temper- 
ature T B ~ 2.76 x 10 9 K. The stimulated emission rate R for 
saturated masers in this transition is given bv lKemball et all 



Figure 4. Enlargement of one feature from the NE of the maser 
shell imaged at <fi = 0.818, from Fig. [2] The vectors show the orien- 
tation of the EVPA and 2 mas in length represents P 1 Jy beam . 
The inner and outer magenta circles have diameters correspond- 
ing to dp and L, respectively. 



(2009): 



R = 23 



2 x 10 10 K 



dQ, 
10- 2 sr 



(6) 



giving R ~ 15 s . 

The Zeeman splitting rate for SiO is given by |^ ~ 
200B s" 1 (|Plambeck et al.ll2003f ). Using the value of B esti- 



mated above, gfl ~ 900 s~ . For this work, the Zeeman split- 
ting ratio is greater than the stimulated emission rate (gQ ^> 
R). This result excl udes intensity-dependent non-Zeeman 
circular polarization (|Nedoluha fc Watsonlll994T ) . The radia- 
tive decay rate V for t he v — 1, J — 1 — transition is < 1 s _1 
l|Kemball et al.ll2009l) . Our data are in the li mit g S7 3> R 3> 
r, which indicates that the lGoldreich et al.l (|l973T ) solutions 
for the linear polarization are applicable. A field of this mag- 
nitude (~ 1 G) also effectively rules out anisotropic radiative 
pumping models on the grounds that the Zeeman preces- 
sion rate is sufficient to mix the magnetic sublevel popula- 
tio ns efficiently. A similar conclusion has also been reached 



by Asensio Ramos. Landi Deel'Innocenti. fc Truiillo Buenol 
l|2005h . 

The ratio of jr ind icates the sat uration level of the 
maser. Equation 8.6.2 of lElitzurl ( | 19921 ) shows that, for our 
SiO transition, the masers are saturated when the bright- 
ness temperature Tb > 0.5 x 10 9 K, confirming that most of 
the R Cas SiO masers are saturated. We therefore interpret 
our results using the magnetic field model for the origin of 
linear polarization. 

Faraday rotation is proportional to the square of the 
wavelength (A 2 ), and so it is smaller at higher frequencies. 
However, Faraday rotation is also proportional to the elec- 
tron density, which is higher at the region close to the star, 
than at higher distances where the longer-wavelength masers 
originate. Previously, it has been presumed that Faraday ro- 
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tation is negligible for wavel ength 7 mm as in our case. W e 
used the relation below from iGarica-Barreto et all |l988h : 

where n e is the electron density and L is the path 
length in the maser region (maser shell thickness). From 
iReid fc Mentenl (jl997h we found that the electron density in 
the SiO maser region ~ 1500 cm -3 , hence the Faraday rota- 
tion is about 16° for L ~ 2x 10 13 cm for a magnetic field of 
B ~ 750 mG. T his estimate is in fac t slightly smaller than 
the estimates by iHerpin et al.l (|2006T ). so Faraday rotation 
could be somewhat higher. It is still within the uncertain- 
ties of our evpa measurements. 



4.2 Linear polarization approaching or exceeding 
total intensity 

We found that in some isolated features, the percentage 
of linear polarization is greater than 100% at three epochs 
cj, = 0.744, <f) = 1.158 and = 1.305,when it reaches 445%. 
This is likely due to resolving out of the total intensity 
Stokes / even on the shortest VLB A baselines. Emission 
on scales greater than ~ 5 mas cannot be imaged by the 
VLB A. Compari s on of the auto- and cross-correlation spec- 
tra (|Assaf et all (|201ll ). Figs Bl-3) shows that 10-90% of 
the total intensity emission of R Cas is resolved out. Fig [6] 
compares the flux density as a function of baseline length 
for Stokes I, Q and U for the spectral range including the 
feature with m< >100% at 0=1.305. Other features in the 
same channels have lower m^ but the average visibility am- 
plitudes still show that the total intensity emission rises on 
large scales (shortest baselines) to a much greater extent 
than polarized emission. The continuum calibration source 
0359+509 does not show any such effect. 

More intense polarized than total intensity has been 
observed before, in low-frequency continuum observations 
using WSRT due to Faraday rotatio n creating smaller-scal e 
structure in the polarized emission (|Haverkorn et alj|2003f ). 
In R Cas, this could be due either to fluctuations in the mag- 
netic field, or to inhomogeneities in the ionisation fraction, 
if these are on smaller scales than turbulence in the neutral 
medium. 

We examined the relationship between auto- and cross- 
correlation spectra for other epochs of high fractional polar- 
ization, for example me ~ 70 — 80% around 27 km s" 1 at 
<t> = 1.158 and 1.305 (Fig. [H? ]) . I Assaf et all (|201lh . figs B2-3 
shows that about 50% of the total intensity flux is recovered 
in this velocity range at <j> = 1.027 — 1.432. This suggests 
that the 'true' fractional polarization could be no greater 
than 30-40% (including uncertainties). A more exact com- 
parison is not possible due to lack of spatial resolution in 
auto-correlation data. 



4.3 90° change in polarization angle 

The linear polarization image of the north-east component 
in Fig shows an abrupt transition in evpa (Electric Vector 
Position Angle) of approximately ^ near Vlsr ~ 26.78 km 
s _1 . At this point, the linearly polarized intensity is near its 



Epoch 


<b 
v 


d£l sr 


T B K 


TSXp 


BD62A 


0.158 


0.066 


1.75x10 s 


22.15 


BD62B 


0.241 


0.044 


0.96X10 9 


24.36 


BD62C 


0.310 


0.039 


0.90 xlO 9 


35.65 


BD62D 


0.390 


0.041 


0.97 xlO 9 


30.10 


BD62E 


0.452 


0.039 


1.18 xlO 9 


30.66 


BD62F 


0.521 


0.044 


1.25 xlO 9 


19.40 


BD62G 


0.595 


0.028 


0.38 xlO 9 


44.97 


BD62H 


0.675 


0.030 


0.75 xlO 9 


52.71 


BD62I 


0.744 


0.048 


1.78 xlO 9 


57.91 


BD69A 


0.818 


0.090 


1.14 xlO 9 


36.82 


BD69B 


0.873 


0.030 


0.13 xlO 9 


11.13 


BD69C 


0.956 


0.083 


12.15 xlO 9 


16.47 


BD69D 


1.027 


0.041 


1.27 xlO 9 


15.74 


BD69E 


1.091 


0.022 


5.13 xlO 9 


21.02 


-oDo9r 


1.158 


0.043 


n now i n9 
U.ZZ X 1U 


24.24 


BD69G 


1.234 


0.023 


2.38 XlO 9 


16.04 


BD69H 


1.305 


0.029 


1.02 xlO 9 


23.88 


BD69I 


1.374 


0.052 


3.89 xlO 9 


21.14 


BD69J 


1.432 


0.040 


1.08 xlO 9 


16.84 


BD69K 


1.501 


0.051 


3.57 xlO 9 


16.09 


BD69L 


1.579 


0.097 


1.19 xlO 9 


20.44 


BD69M 


1.662 


0.062 


0.45 xlO 9 


23.19 


BD69N 


1.783 


0.032 


19.85 xlO 9 


25.84 



Table 2. Stellar phase <f>, error- weighted maser beaming angle 
dQ, error- weighted maser brightness temperature Tb and mean 
fractional linear polarization m^. 



minimum. The orientation of the polarization vectors is tan- 
gential to the projected ring in the part of the feature closer 
to the centre of expansion, whilst they are radial further out. 
Fig.[8]shows the fractional linear polarization for this feature 
as a function of distance from the centre of expansion. Note 
the expanded velocity scale. We ordered the components by 
increasing radial distance from the centre of expansion r and 
found the error-weighted mean linearly polarized and total 
flux densities in ten bins in order to calculate the fractional 
polarization, plotted as a function of r. Each point is labelled 
with the mean evpa and symbol size is proportional to the 
polarized intensity. The emission at r < 25 mas has a mean 
evpa ~ —55° and the emission at r > 25.5 mas has a mean 
evpa ~ +30°. The fractional polarization passes through a 
minimum (~ 10 percent) at r ~ 25 — 25.5 mas, where the 
mean evpa passes through the transitional angles of -74° to 
-90°. The linearly polarized intensity is lower in this range, 
compared with most emission at more extre me evpa values. 

This is consistent with the prediction of lGoldreich et all 

l| 19731 ) when the angle 9f between the magnetic field and the 
line of sight changes from < 55° to > 55°, if the magnetic 
field is radial with respect to t he star. A simila r phen omenon 
in TX Ca m was described bv | Kemball et alj f|201 lh and in 
W43A by IVlemmings fc Diamond! i|2006h . 

4.4 Time variability of polarization 

Fig.|5]shows the relationship between linear polarization and 
total intensity in the whole maser shell and in the inner ring 
(radius 25 mas) . There is an anticorrelation between the to- 
tal intensity peak and the fractional polarization during the 
first optical cycle. The percentage flux density is similar or 
slightly higher in the whole shell compared with the inner 
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45 90 135 45 90 135 45 90 135 45 90 135 
Polarization angle - shell position angle (degrees) 

Figure 5. Histogram of the distribution of polarized emission 
with respect to the deviation of the polarization angle from the 
radial direction (evpa— 8). The line thickness is proportional to 
the logarithm of total polarized emission intensity (at all angles). 
The maxima occur at tf> = 0.390, reaching 160 Jy in the inner shell 
and 350 Jy in the whole shell. The minima occur at <j> = 1.432, 
falling to 4 Jy in the inner shell and 7 Jy in the whole shell. 
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Figure 6. Stokes /, Q and U flux density as a function of baseline 
length for channels averaged from 26.3 - 27.4 km s — 1 , <j> = 1.579. 



shell up to (f>— 0.595. The highest fractional polarization, at 
(f>= 0.744, comes mainly from the inner shell. The fractional 
polarization is also higher in the inner shell from 4>— 0.873 
to (f>= 1.234. At later epochs, the whole shell has higher po- 
larization but at epoch 1.432 and later, most of the emission 
comes from the inner shell. 

We investigated the variation of mean linear polariza- 
tion as a function of stellar pulsation phase in Fig. 1101 We 
calculated these values by averaging over all components 
in each epoch. The error bars represent the standard devia- 
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Figure 7. The NE maser feature seen at (j> = 0.818 in Fig.[2]and 
enlarged in Fig. [4] The mapped emission was summed over Vlsr 
25.67 - 27.43 km.s -1 . The contours are at (-1,1,2,4,8) X 1.7 Jy 
beam -1 . The vectors show the EVPA, plotted every two pixels (0.1 
mas), length 1 mas = 0.8 Jy beam -1 linearly polarized intensity. 




26.5 26.0 25.5 25.0 24.5 24.0 23.5 
Radial distance from centre of expansion (mas) 

Figure 8. The feature shown in Fig. [7] The error- weighted mean 
percentage polarization is plotted for ten increments in radial dis- 
tance from the centre of expansion (r). Symbol size is proportional 
to the mean polarized intensity in each bin, from a minimum of 3 
Jy around r = 25.5 mas, to 11 Jy at r = 23.8 and 26.2 mas. The 
symbols are labelled with the error-weighted mean polarization 
angle. Note the expanded velocity scale. 



tion of average fractional linear polarization. Over the whole 
shell, the mean linear polarization rises from <f> =0.158 until 
4>= 0.744 when it reaches the maximum value. The mean 
fractional polarization then drops abruptly from the maxi- 
mum to reach its minimum at (j>— 0.873. After that, it tends 
to increase again till cf>= 1.783. 

Total intensity, polarized intensity and fractional polar- 
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Figure 9. Distribution of polarized intensity with respect to 
epoch in the inner ring and whole shell 




Figure 10. The relationship between the fractional linear polar- 
ization and the stellar phase. The error bars represent the stan- 
dard deviation. 



ization for the inner ring (25 mas radius), and for the whole 
shell, were plotted with Hanning smoothing over 5 channels, 
giving ~1 km s _1 effective resolution. 

• Fig [EU (Epochs 0.158-0.452), The ratio of fractional 
polarization in the inner shell to the whole shell barely ex- 
ceeds 1. The overall flux density is high at these epochs and 
the fractional polarization is mostly < 50%, higher at red- 
shifted velocities where the emission is almost entirely from 
the inner shell. However, at blue-shifted velocities, where 
there is more extended emission, this is usually more polar- 
ized. At epoch 0.310 and later, the fractional polarization in 
the inner shell is slightly higher than in the whole shell, at 
a few central velocities. 

• Fig [H (Epochs 0.521 - 0.818), The overall flux den- 



sity decreases at all epochs except 0.818. The fractional po- 
larization is quite variable, and at epoch 0.744 the fractional 
polarization from the inner shell exceeds 100%. The ratio of 
fractional polarization in the inner shell to that in the whole 
shell also exceeds 1, and although the uncertainties are high 
the excess is significant, mostly at intermediate velocities. 

• Fig [U (Epochs 0.873 - 1.305), The trend continues, 
with higher fractional polarization from the inner shell at a 
range of velocities. The total intensity from the whole shell 
increases, but with fractional polarization mostly < 30%, 
whereas from the inner shell it exceeds 50% at most epochs 
at some velocities. 

• Fig [TJ] (Epochs 1.374 - 1.783), The total intensity 
decreases, especially at red-shifted velocities, and the emis- 
sion is mostly from the inner shell, so the ratio of polarized 
intensities is close to 1, although it exceeds 1 at epochs 1.374 
and 1.783, when there is more emission from the whole shell 
but it is on average less polarized than from the inner shell. 

Fig. [9] shows that the fractional polarization tends to 
be greater when the total intensity is weaker. Figs. [TT] - [T3] 
show that this is often due to one or two narrow spectral re- 
gions, usually near the centre of the spectra, which have high 
fractional polarization. These are associated with the inner 
shell, see <j> 0.675, 0.744 and 0.956 in the first cycle. Similar 
behaviour is seen in the second cycle at several epochs. 

Thus, there is a distinct tendency for higher fractional 
polarized intensity to be associated with weaker total inten- 
sity and with the inner shell. The emission at more extreme 
velocities tends to come from the inner shell, as projected 
on the sky. This is not surprising since, for a spherical shell 
with a radial velocity field, we would expect to see the most 
red- and blue-shifted emission close to the line of sight to the 
star. Not e that red-shifted e mission within the innermost 25 
mas (R-t, 'Weigel t et ali 2 000l ) must be on the near side of the 
star, seen in infall. 

There are no obviously cyclic repeating trends in the 
second cycle compared with the first. The first half of the 
first cycle shows lower fractional polarization in the inner 
shell for most spectral features whilst the inner shell domi- 
nates the polarization in the second half of the cycle. In the 
second cycle, most of the polarization comes from the inner 
shell but there is no clear dependence on phase. 

What does this imply? If a magnetic field is responsible 
for polarization, this would be expected to be stronger nearer 
the star, and might also be enhanced by shock compression. 
Since the shock crossing time would not necessarily be an 
exact multiple of the stellar period, this could produce the 
out-of-phase enhancement. If, on the other hand, the mecha- 
nism were radiative, that could also be strongest in the inner 
shell but would be expected to follow the stellar phase, so 
this is less likely. 



4.5 Alfven speed 

If we use the estimate of B based on equating the thermal 
energy density with the energy density in the magnetic field, 
we expect the Alfven speed va and sound speed v s to be 
similar. 

We calculated the sound speed by using the relation: 



jkT 



(8) 
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Figure 11. Ratio of fractional linear polarization in the inner ring, to the whole shell for the phases as labelled at epochs 1-5. 



o 



O 



3.0 

2.5 

2.0 

1.5 

1.0 

0.5 

0.0 
3.0 

2.5 

2.0 

1.5 

1.0 

0.5 

0.0 



0.521 




0.595 



illlllllll 



0.675 




+ + + + " ^ + 



H h 



0.744 



0.818 



0.873 



A, 



] \\ I l+ l 



22 24 26 28 30 32 



22 24 26 28 30 32 22 24 26 28 30 32 
Vlsr (km/s) 



Figure 12. Same as Fig llll but for epochs 6—11 
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Figure 13. Same as Fig llll but for epochs 12-17 
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Figure 14. same as Fig llll but for epochs 18-23 
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where 7 is the adiabatic constant, k is the Boltzmann con- 
stant and T is the effective temperature in the masing region. 
We found the sound speed v s = 2.77 km s _1 . The total ran- 
dom energy in the gas is likely to be slightly higher due to 
other turbulence. 

The Alfven speed va is given by: 



where p is the total mass density for well-coupled fluid and 
it is given by 

p = y~\n s m a (10) 

where fj, is the perme a bility of vacuum (4tt x 1(T 7 ) H m 
From iReid fc Mentenl (|l997l ) we found the total density is 
0.4 x 10 n cm~ 3 at the midpoint in the maser shell (7.3 x 10 13 
cm). We found va =3.5 km s" 1 using B = 7 50 mG, or higher 
if the values of B from iHerpin et al.l (J2006) are considered. 
Thus, in the SiO maser clumps, va and v s are similar. In a 
more diffuse, highly-ionised region, with lower density and 
poorer coupling, va could be considerably higher. 



5 CONCLUSIONS 

The SiO maser emission is significantly linearly polarized. 
Most of the linear polarization vectors are either perpen- 
dicular or tangential to the projected shell, but other an- 
gles can be seen. During the first cycle, up to phase 1.158, 
two-thirds of epochs have the majority of the emission with 
polarization angles parallel to radial direction, i.e. parallel 
to the direction of outflow of the first cycle. Emission with 
perpendicular position angles is commonest in the other five 
epochs. Emission with intermediate polarization position an- 
gles dominates the remaining epochs, but the intensity tends 
to be lower leading to greater uncertainties. 

The bi modal distribut i on of evpa is consistent with that 
reported by ICotton et.al.1 (I2006T ) in VLBA observations of 
SiO masers in Mira variable stars at 7 mm. They reported 
that for some features the evpa are tangential while for 
the others are radial. The linear polarization evpa is pre- 
dominantly tangential to the projecte d intensity shell in th e 
VLBA observation toward TX Cam (|Kemball et all [2009). 
This pattern has been confirme d in further observati ons to- 
ward jrJCCainand IRC+ 10011 (jPesmurs et.al.ll2000l '). How- 
ever, ICotton et.al.1 (|2004r i found there was no prevalent pat- 
tern of evpa in any observed stars. 

Our data are in the limit gO R > T, which indi- 
cates that the iGoldreich et~aT1 ll973T ) (GKK) solutions for 
the linear polarization are applicable, i.e. the standard inter- 
pretation for the radially and tangentially evpa orientations 
relative to projected magnetic field is applicable. This re sult 
is consistent with that reported bv lKemball et af] (|2009h to- 
ward TX Cam. The high brightness temperatures observed 
indicate that the R Cas masers are predominantly saturated 
so we consider the GKK model for saturated emission. We 
used the GKK model to estimate the depth, I, above the 
midplane of the shell, of the feature showing an abrupt 90° 
polarization angle change, see Fig 1151 At the point where 
the (dashed) line of sight crosses the vertical magnetic field 
line, in the plane of the sky, the angle between the mag- 
netic field direction and the line of sight, Of, is 90° and 0° ^ 




Figure 15. Cartoon showing an observation of a maser shell with 
a radial magnetic field. The angle between the magnetic field and 
the line of sight 8p is < 55° in the hatched area, producing EVPA 
parallel to the magnetic field; in the solidly shaded area the angle 
is > 55° and the EVPA is perpendicular. At the transition point 
the maximum line of sight through the shell is 21. 



evpa ^ 35° with respect to an axis perpendicular to the lo- 
cal magnetic field direction. When 8p reaches 55° (the van 
Vleck angle), the evpa becomes parallel to the direction of 
the magnetic field as projected against the plane of the sky. 
At this point, Eq. [11] gives I equivalent to ~ 15 mas or 2.6 
AU. 



cos55=- . (11) 

-Rshcll 

If the magnetic field is radial, then the preponderance of 
evpa parallel to the radial direction (during the first cycle, 
with the best quality data, Section T3. 2 [) suggests that a sub- 
stantial fraction of the masers which we detect come from 
a region at least ±2.6 AU deep with respect to the plane of 
the sky containing the star. This suggests a shell thickness of 
^ 4.6 mas or 0.8 AU, which is indeed similar to or less than 
the values of dr given in Table 1. The shell thickness varies 
from about 1/5 to 2/3 of the shell radius, error- weighted 
mean 1/3. The next most common orientation of evpa is 
tangential, suggesting that it is close to perpendicular to 
the magnetic field, when emanating from masers closer to 
the plane of the sky. The magnetic field lines must close at 
larger distances or orientations. 

The percentage of the linear polarization is mi ~ 11— >• 
58% averaged over features for each epoch. The highest frac- 
tional polarization, from the inner shell, is seen in some fea- 
tures in epochs where the total intensity is lower. There is a 
tendency at all epochs for the fractional polarization in the 
inner shell to be greater at the line peaks. Figs. Ill I to 1141 do 
not show any other relationships between fractional polariza- 
tion and velocity. Neither is there any discernible correlation 
between 9f and velocity. 

The changes in the fractional linear polarization 
strength are probably due to fluctuations in the maser am- 
plification process rather than the magnetic field strength. 
The changes in polarization vector orientation are likely 
to be due to small variations in the magnetic field direc- 
tion which is possibly caused by local turbulence (see Sec- 
tion 14. ip . Our estimates of the magnetic field strength in 
the maser region make pumping models based on radiation- 
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driven anisotropies in the magnetic sublevel populations 
highly unlikely. 

A few isolated features with fractional polarization ap- 
parently exceeding 100% can be explained if there is struc- 
ture in the ionised fraction or the magnetic field on scales 
<5 mas but the total intensity emission is smooth on larger 
scales and a greater proportion is resolved out by the VLBA. 
Since the Alfven speed is unlikely to exceed the sound speed 
significantly, and disturbances at either velocity would take 
> 1 yr to cross a 5-mas clump, intrinsic variability is un- 
likely. However, Faraday rotation is not negligible at 7 mm 
wavelength. We found that it is ~ 15° for masers propagat- 
ing through the depth of the SiO shell, using our estimate 
of the magnetic field strength B — 0.725 G based on the en- 
ergy balance, or slightly higher if considering the values of 



B measured bv lHerpin et al 



(20065 • K there are small-scales 
inhomogeneities in the ionised fraction of the stellar wind, 
this will affect the propagation of the polarized emission. 
The small-scale structure seems to be imposed on the maser 
polarization between the emitting material and the observer 
but the effect must occur close to the star where the mag- 
netic field is strong enough. This mechanism is also likely to 
affect other maser features with high observed fractional lin- 
ear polarization, such that although the measured me=ll~ 
58 percent, the higher values might be approximately halved 
if all the emission was measured by the interferometer. We 
are, therefore, probably within the limit of the GKK model, 
wherein the mean fractional linear polarization can be up to 
33 percent. 
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